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Abstract 

The ground state and low T behavior of two-dimensional spin systems with discrete binary 
couplings are subtle but can be analyzed using exact computations of finite volume partition 
functions. We first apply this approach to Villain's fully frustrated model, unveiling an unexpected 
finite size scaling law. Then we show that the introduction of even a small amount of disorder on 
the plaquettes dramatically changes the scaling laws associated with the T = critical point. 
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Two-dimensional frustrated systems typically have a high ground-state degeneracy: this 
can lead to unusual physical properties In the absence of disorder, the zero-temperature 
spin states can often be mapped [sj to Baxter type or solid on solid models, leading to a 
relatively good understanding of the low temperature regime. Unfortunately, the incorpo- 
ration of disorder renders analytic approaches powerless and one must resort to numerical 
treatment. In the limit of very strong disorder, one expects to have a spin glass, at least at 
zero temperature. It has recently been shown that the critical behavior of two-dimensional 
spin glasses with binary quenched random couplings is rather subtle, being very different 
from the one suggested from the naive low temperature expansion 0, U, 0, |y] . In this work, 
we focus on the change of critical behavior of the 2d fully frustrated model (FFM) when one 
adds disorder. In the FFM each elementary plaquette of the lattice is frustrated, and the 
couplings between the spins are all of the same magnitude. At T = (the model's critical 
point), the FFM exhibits a power-law decay of spin-spin correlation functions. Furthermore, 
for T > 0, the correlation length diverges exponentially in 1/T. Suppose now we allow some 
plaquettes to be unfrustrated; does this reduction in the frustration leave the critical prop- 
erties invariant, or is it instead a strong perturbation? Using exact partition functions on 
finite lattices, we shall show that even small amounts of disorder dramatically change the 
thermodynamic singularities of the model. 

Villain's Fully Frustrated model - We consider a two dimensional square lattice with 
Ising spins on the sites and binary couplings Jij on the bonds; without loss of generality, we 
set \Ja\ = J = 1. The Hamiltonian is 



where the sum runs over all pairs of nearest neighbor sites and boundary conditions are 
periodic. We shall consider both L x L lattices and the strip geometry (L x oo). In the 
FFM the product of the four couplings on each elementary plaquette is —1: there are many 
different ways to do this and most of them are gauge equivalent. More precisely, for a 
fully frustrated lattice with periodic boundary conditions, there are four equivalence classes, 
corresponding to whether the product of the J along a loop winding around a direction of 
the lattice is +1 or —1. Following Villain we choose the periodic implementation where 
all Js are set to 1 except on the vertical bonds where lines of couplings alternate between 
+ 1 and —1 values. Villain applied the transfer matrix formalism to extract the free energy 




(1) 



2 



density in the infinite volume limit: 

-PfooiP) = ln(2 cosh(/3 J)) + JL (2) 

/ d/i / rfA; ln[(l + 2 2 ) 2 -22 2 (cos2/i + cos2fc)] . 
Jo Jo 

Here z = tanh(/3J), /3 = T _1 is the inverse temperature, and is the free energy per site 
in the infinite volume limit. 

The free energy turns out to be analytic everywhere except at T = 0. The T = entropy 
density of the system (normalized to the number of sites), sq, is finite: so = C/n, where C 
is the Catalan number^, [|. The ground state energy density e is — J (one quarter of the 
links are unsatisfied), and the low T expansion of Eq. ((2)) to leading order is 

Pf oo {P)~/3e -s Q + c 1 l3e- 4f>J . (3) 

This makes clear the non-analyticity of the free energy at the critical point T = 0: a well 
behaved low temperature expansion of flf^ would give a series in y = exp(— 4/?J) while in 
expression (J3J) a factor j3 (i.e., a factor ln(y)) multiplies the leading contribution y = e~ 4/3J . 

Taking the second derivative with respect to T gives the leading behavior of the specific 
heat density: 

cy(/3)~16ci J 2 (3 3 e~ 4/3J . (4) 

In the present case of the FFM, Eq. (J2J) tells us everything about the temperature scaling 
of thermodynamic quantities. But it is also possible to obtain the free energy in a finite 
volume: we shall use this to extract both the correlation length and the finite size scaling. 

The FFM in a finite volume or on a strip - Computing the partition function of the 
FFM in a finite volume is not difficult. For a finite lattice with periodic boundary conditions, 
the partition function associated to the Hamiltonian of Eq. Q can be expressed as the sum 

nnn 

|1(1I11LI12|. each multiplied by a plus or minus sign: because of these signs, 



of four Pfaffians 



there are large cancellations when combining the contributions of the four Pfaffians. Thus 
it is necessary to compute each Pfaffian to very high precision; we have done so with the 
"Mathematica" program that allows for arbitrary precision computations. Basically, in this 
approach, the double integral of Eq. (J2J) becomes a sum over a discrete set of momenta, 
where the support of this set is different for the four Pfaffians. We have also considered 
a strip geometry, where the vertical direction (where couplings are ±1) is infinite and the 
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other has a width L. In this case we have to compute mixed sums and integrals from which 
we obtain the free-energy density of the FFM in the strip geometry (typically we get 15 
or more significant digits). We now focus on how the thermodynamic limit is reached with 
increasing L. 

The correlation length - The diverging correlation length of a physical system ap- 
proaching criticality can be defined in different ways: for example one can use the exponen- 
tial decay of spin-spin correlations or the exponential convergence of the free energy density 
with increasing lattice size. Since in our approach we do not have access to the spin-spin cor- 
relation functions, we use the second definition. Consider the strip geometry associated with 
a lattice of infinite length in one direction and of width L in the other (with periodic bound- 
ary conditions): the correlation length £ can be defined via the relation (/x — f^) ~ e~ L ^, 
where the exponential can have a prefactor that depends smoothly on L, for example via a 
power law. 

In Fig. (JTJ) we show ln(/oo — /l) as a function of L (points) and the best fits to the 
form a(T) — m(T)L — c(T) In (L) (continuous lines), for five temperatures (T G [0.4,0.6]). 
m(T) = 1/£(T) is the inverse correlation length. The quality of the fits is excellent, and we 
find that c(T) depends only weakly on T and is close to 1.5. Because of that, it is appropriate 
to parametrize the L dependence via the correlation length £ through the relation 

e -L/( 

h-fco = A(T) -jjj- , (5) 

where A(T) is a smooth function of T. From this we can extract £ when L is large; since our 
numerical values of fi — foo loose precision beyond L = 200, and getting higher precision 
would have demanded a very large computational effort, we have been able to extract £ by 
curve fitting only for T > 0.4. In the inset of Fig. (0) we plot ln(£) as a function of 1/T, and 
find a straight line with slope very close to 2. Similarly, the prefactor of the exponential is 
very close to | and we conjecture that this is in fact the exact value (note that this prefactor 
has not been estimated previously). We thus conclude that 

£(T) ~ e^/2 . (6) 

This scaling form applies as T — > 0, but it holds to good accuracy even when £ is not so 
large: for example at T = 1 where £ ~ 4 the value given by © is only 10% off from the 
actual measured £. In the standard hyperscaling framework, the singular part of /Sf^ is 
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FIG. 1: ln(/oo — f£) as a function of L (points) and the best fits (continuous lines). In the inset, 
ln(£) as a function of 1/T. 

given by for a d-dimensional model (up to constants and possible logarithmic terms in 
£). Using Eq. (0), we then expect the singular part of to go as exp(— 4/3 J), which is what 
was found from Eq. (jH)): in the FFM, hyperscaling holds. 

Finite size scaling - We have discussed the 1 < ( < L region. It is also of interest 
to consider the regime where L <C £: only by controlling this region we can reach a full 
understanding of the finite size scaling of the system. 

As a start, we notice that good data collapse is obtained as T — ► when we multiply 
(/l — foo) by the factor PL 2 , i.e., we find 



(3L 2 {f L - U) ~ W{L/i) 



(7) 



where W(L/£) is an adimensional function of the ratio L/£. When x = L/£ — > oo we recover 
the previous analysis and W(x) ~ t/xexjp(—x). In the opposite limit, x = L/£ — ► 0, the 
function W goes to a constant Using the relations 0, © and (0), we find 

S -ARl W 



h 



e ~ 



(3 ' ~ x ~ ' (3L 2 ' 

Dividing and multiplying the last term by £ 2 and substituting 2(3 by ln(£) (cf. Eq. ©) we 
obtain 

h * e - | + Cl e" 4 ^ f 1 + W / ln(0) (8) 

where W is a scaling function simply related to W. The important point is that in Eq. (jHJ 
the adimensional scaling function W(L/£) is divided by a factor /3 ~ hi(£), so one has an 
anomalous finite size scaling law. 
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The function W satisfies 

~ x~ 2 when x — > , 

W(x) = < 

I constant when x — > oo . 

We show all of our data for /3((f L (P) - e + f )e 4/3J / c i - 1) versus L/f in Fig. ©• Here 
10 < L < 180 and < T < 1.0. C\ has been determined with a best fit to /oo and has the 
value ci = —1.273. The data collapse is excellent and since it involves a correction to /l, 
one can conclude that finite size effects are under very good control. 

To complete our study of the FFM without disorder, we consider finally the low temper- 
ature expansion of 

W) * eo - ~ -^exp(-4/3J) , (9) 
p p-L g 

where go an d g\ are the respectively the degeneracy of the ground state and of the first 
excited state. We determine go and g%, finding that to very good accuracy 

gi /g ~AL 2 + BL 2 \n(L) , (10) 

with A = —0.44 and B = 0.63. It is possible to show that the scalings (fTU|) and (jHJ) are 
mutually compatible. 

Adding disorder / diluting frustration in the FFM - At this point the large L and 
£ behavior of the FFM is well understood. Now we move on to see what happens when 
frustration is partly removed. We do this by unfrustrating a small fraction of the plaquettes, 
choosing these at random. The set of couplings Jij entering the Hamiltonian (0) is now such 
that on a fraction pi of the plaquettes their product is equal to 1. We shall refer refer to 
this as the Plaquette Disorder (PD) ensemble. Does this modification to the Hamiltonian 
change the scaling laws of the FFM? Interestingly, the change is in fact dramatic, as we now 
reveal. 

The computational tool - The presence of plaquette quenched disorder breaks the trans- 
lation invariance of the Hamiltonian; because of this, the transfer matrix cannot be diago- 
nalized by going to Fourier space. Instead, we rely on the explicit computation of Pfafnans; 
this can be done for any set of J^. A very effective approach for computing such Pfafnans, 
based on modular arithmetic, has been proposed and implemented in Q|. One evaluates 
the partition function in a (large) finite volume via its low temperature series: 

Zj W = e 2L2/3J Pj (e- 2 ^) , (11) 
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FIG. 2: P((Jl(P) - e + f )e 4 P J /ci - 1) versus L/£. Here d = -1.273. Inset: -T\n(T 2 c v ) vs. T 
showing the convergence to A = 4 (cf. section with disorder). 



where Pj(x) is a polynomial whose integer coefficients are the number of spin configurations 
of a given energy. The algorithm determines these integers exactly, allowing one to 
analyze the system even for very low temperatures. 

In our implementation P. Il2j] the CPU time to compute Zj grows approximately as L 5,5 . 
We have mainly studied the model with an unfrustrated fraction of plaquettes pi = 1/8, by 
determining Zj on lattices with sizes ranging from L = 24 (2000 samples) up to L = 56 (200 
samples) . We have also analyzed the case of p± = 1/4 on lattices with L = 32 (600 samples) 
up to L = 48 (100 samples). We computed sample averages of different physical quantities 
like the free energy and the specific heat, and we analyzed T = quantities like the number 
of ground states and of low-lying excited states. 

Low temperature scaling of cy - In the following we mainly analyze the specific heat, 
Cy, that we compute from the fluctuations of the internal energy 

In the case without disorder, the low T thermodynamic limit behavior (where V diverges 
at fixed low T) of cy is given in Eq. (J3J). On the contrary in the finite size limited, T — ► 
limit (where T — > at fixed volume), we have that 

16/3 2 J 2 <7i _ 4/3J 



cv 



(H)} 2 ) 



(12) 



L 2 go 

Thus, the scaling in the thermodynamic limit differs from the finite size limited, finite 
volume low T scaling, but only in the power of the f3 prefactor, and not in the argument 
of the exponential. (This power is related to the logarithmic corrections in the L scaling of 
the ratio gi/go of Eq. (|10))). Disorder completely changes this picture as we shall now make 
clear. 
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Let us parametrize the scaling behavior by the argument of the exponential function and 
by a power for the first subleading correction: cy ~ [3 P e~ A/3J . We have seen in the FFM 
that A = 4 and that P = 2 in the naive, finite size limited scaling while P = 3 in the scaling 
regime of the thermodynamic limit. In Fig. (j3J) we plot —Tin (T p cy) versus T for the PD 
model (with pi = 1/8). We take here P = 2, and try to determine A. We proceed in this way 
since the value of A is not very sensitive to the value that we assume for P: our conclusions 
for the value of A would be unchanged by taking P up to 10, and would be strengthened 
by taking P < 2. A is given by the intercept of the envelope's extrapolation to T = 0. We 
can distinguish three regions in Fig. (jHJ). The region at very low T values corresponds to 
the naive scaling with A = 4 as can be seen from the T = intercept; this region shrinks 
to zero with increasing lattice size. There is also the high T region but it is not relevant for 
critical properties. Finally, the most interesting region is in between the other two; there, 
one has the scaling in the thermodynamic limit, obtained from the envelope of the set of 
curves. It is clear from the figure that this scaling is incompatible with A = 4; in fact, the 
estimated value for A decreases with the lattice size (since for higher L we can determine 
the envelope to lower T values), and possibly goes to zero when L — > oo: such a behavior 
would imply an algebraic T — > limit, and not an exponential singularity (it is clear that if 
this happens P will become smaller than zero: our best fits are already compatible with a 
small negative value of P, but since they are not very sensitive to the value of P we cannot 
give a quantitative estimate of this effect). We have repeated this analysis for a dilution of 
1/4 with similar conclusions: again the curves go down to low values of A, possibly 0. As 
a final note and to drive home even more the incredibly strong effects of disorder, one can 
compare the behavior of cy in the pure and the PD models: as shown in the inset of Fig. (0) 
the pure model has both a very clear exponential low temperature limit and small finite size 
effects; on the contrary in the presence of disorder, cy has large finite size effects and the 
exp(— 4/?J) scaling is clearly absent. 

The degeneracy of low energy states - In the insets of Fig. @ we plot gi/goL 2 as a 
function of L for the FFM and \n(gi/ goL 2 ) (that we get as an output of our exact computa- 
tion of Z) as a function of L for the PD: this is because the first ratio grows logarithmically 
in L (see Eq. (HIIJ)), while in the presence of quenched disorder, it seems to grow at least 
exponentially fast with L. The plaquette disorder completely changes the scaling of these 
quantities and so the low temperature expansion for the pure and the disordered models will 
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FIG. 3: —T\n{T 2 cv) versus T in the PD model with a fraction 1/8 of unfrustrated plaquettes. In 
the inset: on the left, \n(gi / g^L 2 ) versus L for the PD model; on the right, gi/goL 2 versus L for 
the FFM. 

be very different. 

Summary and discussion - Two dimensional frustrated models, with or without dis- 
order, are challenging systems, especially when they have a critical point as Villain's fully 
frustrated model does. We first studied in depth his (pure) model: beyond a conjecture for 
the magnitude of the correlation length, we showed that finite size scaling was anomalous. 
We also found that the ratio gi/go of the degeneracies of the lowest energy states grows as 
L 2 with multiplicative logarithmic corrections: as a result, the power law in T multiplying 
the exp(— 4/3 J) scaling of the specific heat cy is modified. 

We then introduced quenched disorder in the form of a small fraction of randomly posi- 
tioned unfrustrated plaquettes to find that qualitatively new phenomena arise. For instance, 
gi/go grows exponentially in L rather than as a power. Following the argument of what oc- 
curs in the pure model, this growth breaks the exp(— 4/3 J) scaling of cy, taking one to a form 
of the type exp(— Af3J) with A rather small if non-zero. These effects are striking and show 
the extreme fragility of the pure system: the universality class of the FFM is completely 
changed when disorder is introduced. 
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